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H ■ Abstract 
03 

The finite XXZ model with boundaries is considered. We use the Matrix Product Ansatz 
(MPA), which was originally developed in the studies on the asymmetric simple exclusion 
process and the quantum antiferromagnetic spin chain. The MPA tells that the eigenstate 
of the Hamiltonian is constructed by the Zamolodchikov-Faddeev algebra (ZF-algebra) and 
fSJ ' the boundary states. We adopt the type I vertex operator of U q {sl2) as the ZF-algebra and 

realize the boundary states in the bosonic U q {sl2) form. The correlation functions are given 
\ by the product of the vertex operators and the bosonic boundary states. We express them 

■ in the integration forms. 

00 1 1 Introduction 

Even for the integrable systems, to calculate the Appoint correlation functions is a very difficult 
problem. In the XXZ model, the correlation functions are found for the infinite Jl|, |2j, ||| and 
the half-infinite [Q, || chains by the approach of the vertex operators || M. Those vertex 
operators are constructed in the bosonic U q (sl2) form and are known to have two types, type 
I and type II. The former means the half-infinite product of the R-matrices, and the latter 
describes the excitation of the model. In the infinite case, it is natural that the quantum affine 
algebra U q {sl2j is adopted, because the model has its symmetry, i.e., the Hamiltonian commutes 
with its generators. However, the half-infinite XXZ model obviously has no such symmetries. 
In analysis of that model, it is important to know the physical meaning of the type I vertex 
operator. The transfer matrix can be defined by the type I vertex operator without physical 
conflicts; the boundary states are realized in terms of U q (sl2) and the correlation functions are 
calculated by using its expression. Since the vertex operators do not indicate the finite lattices 
physically, in the finite chain with boundaries things go wrong in such a way. 

Without regard to calculation of the correlation functions, the analysis of the models with 
boundaries is in general more difficult than that of the models without boundaries. Roughly 
speaking, there are three theories to treat the boundary. The first uses the boundary K-matrix 
in the quantum integrable systems pi 0, 10]. The K-matrix describes the integrable boundary 



conditions and satisfies the reflection equation (see eq. fl2.1|) ). Several solutions are reported 
[11, 12, Dl, 14 1 . The second is the boundary S-matrix method in the integrable field theories 
[15, 1£, 17), IS | . In such field theories, for instance, the conformal field theory, the scattering 
matrix of particles is factorized into a product of two-particle S-matrices, and the boundary 
scattering is described by the boundary S-matrix. Those S-matrices are known to have similar 
properties as the R-matrix and the K-matrix. 
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The third, called the Matrix Product Ansatz (MPA), is different from the previous two 
methods. It was originally developed in the studies on the asymmetric simple exclusion process 
[19, [H], 21] and the quantum antiferromagnetic spin chain [^, [2^, 24]. The asymmetric simple 



exclusion process is a stochastic process where particles diffuse asymmetrically, and are injected 
and ejected at boundaries. Under the MPA, the state on each site is represented by the infinite- 
dimensional matrices determined by the model, and the correlation functions are defined by 
the products of those matrices. Calculating them is formally possible. However it is generally 
difficult to treat such infinite-dimensional expressions. For removing this difficulty, the MPA 
has been recently extended by using the terminology of the integrable systems |25|]. It has 
been proved that the infinite-dimensional matrices correspond to the representations of the 
Zamolodchikov-Faddeev algebra (ZF-algebra). That is an algebra where the generators are 
"scattered" mutually by the R-matrix and at the boundary by the K-matrix. In addition it 
has been shown that a certain eigenstate of the Hamiltonian is given by the product of the ZF- 
algebra and the boundary states. The remarkable point is that the symmetry of the model is not 
required in spite of using the ZF-algebra. This property is the mirror of the original one. The 
ZF-algebra and the boundary states are realized by comparing the result with the known one 
in other theory. By using them, the eigenstates can be computed irrespective of the symmetry 
of the model. It is also possible to find the correlation functions. The MPA thus offers a new 
powerful tool to analyze the model with boundaries. 

The aim of this paper is to calculate the correlation functions of the finite XXZ model with 
boundaries by applying the MPA. We note that the type I vertex operators can be regarded 
as the ZF-algebra, because they satisfy the same commutation relations. We adopt the type 

1 vertex operator as the ZF-algebra, which is justified by the agreement of the result and the 
known one in the half-infinite chain limit. While the vertex operators indicate the half-infinite 
products of the R-matrices, such physical meanings are hidden and only those mathematical 
properties are observed. The reason is that the vertex operators are used only for realization of 
the eigenstate through the ZF-algebra and are not related with the symmetry of the model. The 
detail is discussed in Section |2[ We realize the boundary states in the bosonic U q (sl2) form in 
Section ||. Based on those bosonic realization, we calculate the 2A?-point function that generates 
the correlation functions in Section ||. In Section [| we discuss our results and note some open 
problems. The technical details are summarized in Appendices |A], |B| and |C]. 

2 Matrix Product Ansatz 

2.1 Formulation of Matrix Product Ansatz 

In this section we formulate the MPA by use of the theory of the integrable systems [25[| . We 
consider the model with boundaries. Let R{Q be the R-matrix and K(Q) be the boundary 
K-matrix. They satisfy the reflection equation 

K 2 ((2)R2i (CiC 2 )^i(Ci) J Ri2(Ci/C 2 ) = ^2i(Ci/C2)^i(Ci) J Ri2(CiC2)^ 2 (C 2 ). (2.1) 

In customary £ is called as the spectral parameter. The lower labels attached to operators 
indicate the lattice sites where they act. The transfer matrix is defined by 



T(Cat,...,Ci)=^(Cat 



^ N-l 

II RiMd<N) I ^v(Gv) I II RnACn/Q) I • (2.2) 



i=l / \ i=l 
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The superscripts + and — of the K-matrices refer to the left (site N) and right (site 1) boundaries. 
The ZF-algebra is a set of operators F(() and F*(Q satisfying the following relations [jD], |lfj, 



17, 



Rl2((l/(2)F 1 (( 1 )®F 2 (( 2 ) 
iMCl/C2K(Cl)®^2*(C 2 ) 

K-(()F(()-W 
K + (()F*(() • V* 
W* ■F*(C -1 )lf~(C) 



F 2 (( 2 )®F 1 (( 1 ), 

F^C 2 )^FUCi), 
V-F(C), 
F(C l )-W, 
F*^ 1 ) ■ V*, 
W* -F*(C). 



(2.3) 
(2.4) 
(2.5) 
(2.6) 
(2.7) 
(2.8) 



Here V and W are the left and right boundary states, and V* and W* are those dual states. 
The MPA claims that a product of the ZF-algebra and the boundary states gives the eigenstate 
of the transfer matrix whose eigenvalue is 1. Namely, the "state" 



|Gv, • • • , Ci) = |V • F N (( N ) ® • ■ • ® Fi(Ci) • W), 
(Ci, . . . , Cjv| = (W* ■ Ff (Ci) ® • • • ® ift(Ov) • V*| 



(2.9) 
(2.10) 



is the eigenstate of the transfer matrix. We make the ansatz that the states (|2.9D and ( |2.10| ) 
have mathematical meanings. The proof is straightforward. Using the relations fl2.3| )-(^lj|), we 
have 



T(Ctv,...,Ci)|C7V,...,Ci) 

/ 



II RiACM k n((n)\C 



N-l, 



• CiXn) 



1=1 



J\ Ri,N{(i(N) \Cn-i, ■ ■ ■ >Cl>CjV ) 



i=l 



^(Cat)IGV--,Ci> 
IOv,...,Ci)- 



Similarly, we can show that 



<Ci,...,0v|T(0v,...,Ci) = (Ci,---,C 



A' 



Observing the relations 



T(l, 



1) 



1, 



(2.11) 

(2.12) 

(2.13) 
(2.14) 

#11,..., 1) =0, (1,...,1|# = 0. (2.15) 

We conventionally call the states (^^) and ( |2. 10|) as the stationary state and the dual stationary 
state. In particular cases the stationary state corresponds to the ground state. For instance, in 
the half-infinite XXZ model, the stationary state is known to be the ground state when \r\ < 1 



a 



oc 



N 



T(Gv,...,Ci) 



=d=i 



we further have 



const, x H, 
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(See eq. ( [2.19| ) for the definition of the parameter r.) The stationary state is one of the 
candidates for the ground state. 

Note that the proof ( |2.11 ) is similar to the derivation of the equations that the correlation 
functions satisfy in the half-infinite XXZ model ||. The difference between them is the definition 
of the boundary states. While in the usual theory the boundary states mean the eigenstate of 
the Hamiltonian, in the MPA they indicate the boundaries themselves. The eigenstate in the 
MPA depends on the ZF-algebra and the boundary states, and the model is defined by the 
R-matrix and the K-matrix. It is important to notice that the boundaries and the Hamiltonian 
are treated independently. 

We indicate that the realization of the ZF-algebra is independent of the number of the sites 
of the model. In the MPA the ZF-algebra is defined by an interaction on the nearest sites and 
the number of the sites of the model is given by the number of the ZF-algebras appeared in the 
stationary state (see (|2.9|)). Then if the ZF-algebra is realized for a particular number of the 
sites, such realization is always correct. 



2.2 Application to Finite XXZ Model with Boundaries 

Let us apply the MPA to the finite XXZ model with boundaries. The Hamiltonian of the model 



is 



7V-1 



i=l 



where erf (a = x, y, z) are the Pauli matrices acting on a site i, 



' 


1 " 




V 

°i = 


' -i " 




°i = 


h- > 


" 


1 





i 


i 


1 

% 





-1 



(2.16) 



(2.17) 



The R-matrix and the K-matrix are defined in Appendix [A|. The derivative of the transfer 
matrix gives the Hamiltonian, 



d 



N 



■T(Gv,...,Ci) 



4q 



<i=1 



(2.18) 



The constants A and h are parameterized as 



A 



q + q- 



l + r± 



Aq 1 — r± 



(2.19) 



We consider the model in the massive regime A < — 1, i.e., — 1 < q < 0. In the MPA, one must 
realize the ZF-algebra and the boundary states satisfying the relations (|2.3| )-( |2~8|) . We adopt 
the type I vertex operator ^(C) of U q (sl2) as the ZF-algebra and assume that the boundary 
states can be described in terms of U q (sl2)- (The vertex operators of U q (sl2) are formulated in 
Appendix [B].) By the definition Q2,9| ) the stationary state is defined by 



|Civ,...,Ci) = VGx |V-$(6v)® 



*(Ci)-w>, 



(2.20) 



with some normalization factor G. By using the dual vertex operators $*(C) ; the dual stationary 
state is similarly defined by 



(Ci,... ) Cat| = VG'x <W*.$*(Ci) 



$*(0v)-V*|. 



(2.21) 
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With these definitions we express the correlation function involving the local operator O by 



(Ci,...,0v|O|0v,...,Ci> 

= G x (W* • **(&) ® • • • ® $*(Gv) • V*|0|V • $(Gv) <8) • • • ® *(Ci) • W). (2.22) 

The definitions ( [2.20D and ( |2.2lD are justified by the verification of the expression (|2T2l ) in the 
limit N — > 00. We assume that the local operator (9 has a 2 x 2 matrix form = J2 m =i(^rn)ee' ■ 
In the limit N —* 00 the boundary states V and V* are ignored. Applying the invertibility of 
the vertex operators (B.13) we find that 



(Ci,...,0v|O|Gv,...,Ci> 

TV 

= g m ~ N G X E E < W * • ^(CO • • • ^L.SCm-lMJCm) 
m=l ei,...,e ra ,e{ n 

^(OmJ^^yVx^i) • • • 3* (Cl) • W), (2-23) 



B.14). We have used the vector expressions of the 



where g is a scaler constant given in eq. 
vertex operators $(C) and <3?*(C) (see eqs. ( B.5| ) and ([B.1CQ ) 

$(C) = E * e (0 ® w e , **(C) = E $ *(0 ® <• (2.24) 

We set G = g N /(W* • id|id • W). If (W* • id| and |id • W) are the same as the eigenstates of the 
half- infinite XXZ model, the correlation function ( 2.23; ) is identical to the already known result 



(Notice that our boundary states correspond to the eigenstates in the usual theories.) In the 
next section we shall compute the boundary states and make sure that this conjecture is indeed 
correct. The correlation functions correspond to the known results in the limit N —* 00. Note 
that the realization of the ZF-algebra is independent of the number of the sites of the model. 
The justification is thus approved. 

Remarks are in order. First, in the present case, U q {sl2) does not play the role of a symmetry, 
because it does not commute with the Hamiltonian ( |2.16| ). We adopt the vertex operators for 
realization of the MPA and disregard the mathematical symmetry of the model. It is not 
requested in the MPA. Second, one might consider that the physical meanings of the vertex 
operators conflict to finiteness of the model. This doubt is reversed by the following reasons. 
In the usual theory the vertex operators constitute the transfer matrix and prescribe the model 
through it, however, in the MPA they realize the ZF-algebra and give the eigenstate but not the 
transfer matrix. Only the R-matrix and K-matrix compose the transfer matrix, and the vertex 
operators do not appear in. The roles of them are difference. Finally, we point out an open 
problem of the MPA. We have adopted the vertex operator as the ZF-algebra only because the 
result must be identical to the known one in a certain limit. For application of our MPA, one 
needs results of other theories to realize the ZF-algebra. There are yet no prescription to derive 
it within a scope of the MPA. 

3 Realization of Boundary States 



In this section we realize the boundary states in terms of Uq^sl?). The relations (|2.5| )-(2.< 
request 

v-$(r 1 )^ + (C) = v-$(0, K-(CMO-w = HC 1 )-^, (3.1) 
k + (CW(0 • v* = $*(C -1 ) • v*, w*-$(C _1 )if (C) =W*-$(C). (3.2) 
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We make the ansatz that they have the forms 

V = e FiV) ® e- A -/ 2 E , W = e F(W) ® e A E ftC ) , (3.3) 

y* = e*" 00 ® e A+ «/ 2 E H« W* = e F( "*' ® e" A E W<°>. (3.4) 

The definitions of the spaces (i = 0, 1) and the operators e A , e a are given in Appendix ^. 
Since the spectral parameters are reversed £ — > £ _1 by the action of the K-matrices, the prior 
parts of the boundary states i* 1 ^, F^), F^ v *^ and i* 1 ^*) must have the quadratic terms of 
bosonic operators. They are expressed by the following forms, 

1 oo 7 (V) oo 









1 oo 

2 ^ 

fe=l 


, (W) oo 

»] " k h 


(3.6) 








-, oo 

= ~Y 

2 ^ 

fe=l 


r (V) oo 


(3.7) 






F (W) : 


1 oo 

- ^ 

z fe=l 


, (W*) oo 

[2k]\k] ak+ h k 




where [re] is the 


q-integer 


• (see eq. 


(|B.2j)). Noting the following formulas 






e A e B 


= e ±[A,[A,B]] e lA,B] e 


B e A , (for [A, [A, B]] E C) 


(3.9) 




e A e B - 


= e B e A e [X,B] e I[[AB],B] ) (for [[A 5J B],i?] E C) 


(3.10) 


and the relation 


ak • 1 = 


1 • a_ fc = 


(fc > 0), d ■ 1 = 1 • d = 0, we find that 








a {V) 


= -oik, 




(3.11) 






(W) 

a k = 


= -at 1 , 


pP = 


(3.12) 






(V) 

a k 


= -at, 


p^ ^ = /5fc+> 


(3.13) 






(W*) 

a k = 


*-l 
-«* » 




(3.14) 


where 














«fc = 


<T 6fc , 


Pk± 


„-7fc/2 

= i^r ( ^ (1 " 9fc)+gfcr ± ) ' 


(3.15) 



k /2 

= g 2fc , = -p^ " Q k ) ~ q~ k r± k ), (3.16) 

and Ok = (1 + (— l) fc )/2. The relations ( |3.11| )-( |3.14| ) are easily proved by the following identity 
(see eq. ( |A.2| ) for the notation (z;p)oo), 

The boundary states W and W* are equal to the eigenstates in ref j|. The MPA has been 
approved. 
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4 Correlation Functions 



To find the correlation functions we compute the following 2iV-point function, 

-^e* r-> e ivi e Wr-->ei ({Cj }i {Ci}) 

= G x (W* • <D* (0 '%*) • • • Kt^iCh) ■ V*|V • ^t L) «N) ■ ■ ■ <^' 0) (Ci) • W). (4.1) 



"1 C JV 



Here t = N mod 2 and G = g N /{W* ■ id • V*|V • id • W). Note that the product of the boundary 
states V* and V can be regarded as an element of either or IIS 1 } by the cancellation of the 
operators e a . First we compute the normal-ordering of the bosonic operators. By the bosonic 



formulas (B^), (|B.3| ) and the expressions of the vertex operators (BJ3) and ( B.7 ), we have the 
following expressions. 

^-^(^•••^(Ci) = n / ^({ca,w)x/({ca,w), (4.2) 

dw 



Kf' 1) (C 1 )---Kt L,t) ^N) = nf,^^({Ci}»W})x f ({Ci}.K})» (4-3) 

beB c b Z7n 



where 



? 4^2/, 2.,, I 



(4.4) 



CW.KH = (- g 3)l[f]^-E^ & n(- g - 1 c;)^ + ^ i n i|or 



Ilbgj? q Hi - g 2 K IWK - OK - g 2 Q 

ni< 6 t 2 (c* 2 - o ni>6K - q 2 c; 2 ) 



(4.5) 



oo 



J({Ct}>{«>a}) = exp I ^a_ fc X fc ({Ci},{w a }) 1 exp I -^a fe yfc({Ci},{w a }) j 

\fe — 1 ' » fe — 1 / 

8 e-^^nH 3 C^ /2 II< S . ( 4 - 6 ) 

i=l aGA 

r({c;},K» = «ap ^^({^.M})! (-E"*y;({{;}. W})) 

AT 

o^i^ n(-</c; 2 ) a/2 n < -a > ( 4 - 7 ) 

i=l beB 



with 



„7fc/2 N k /2 

**({<*}, W) = ^y-ECf-^E^, (4-8) 



-5fc/2 A 7 fc/2 

n({6},W) = ^Ec 2fc -^L>a fc , (4.9) 
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„3fe/2 N k/2 

xmc*hM}) = ^Ecr-^E«**. ( 4 - 10 ) 

-k/2 N k /2 

In the above, sets ^4 and B are denned by 

A = {i\e t = +}, B = {j\e* = -}. (4.12) 

The integration contours C a and in eqs. ( |4.2[ ) and ( [4.3| ) are chosen to encircle the only points 
q 4 d (* < °) an d (? 2 Cj* (i > respectively. By the formulas ( |3.9[ ) and ( |3.10| ), it follows that 

|v •/({&}, W)-w> 

x |e F(V) e F(W) exp I ^ a_ k (X k + a^Yfc) ) exp I ~ E afcYfc ) 

\ k — 1 / \ k — 1 / 

^e-*^^ n(-5 3 c 2 ) 9/2 n v>, ( 4 - 13 ) 

i=l aeA 

(w*-i*({c*},K})-v*| 

N 

*-l 



n(-^ 2 ) i/2 n" 



x <e^ 'e^ J exp f £ a_ fc (X fc * + a^*) j exp ( - £ ^1 

0e t (Ef =1 eJ+D jj ( _ gC *2 )S /2 -q ^-^ (414) 
i=i fees 
Combining these expressions together we obtain 

Pel,...,e* N ;e N ,...,ei{{(jY,{(i}) 

= ^ n n £ L^^cm,{*a})c*({GhM}) 

aeAbeB JCa JC b 

N \-i(Ef =1 ^+i) 

n (-«c; 2 ) n >< w; {c*i, k», (4.15) 



where 



J({< ( }. {»«};{<;},«}) 

= « x e*p {- £ MM ('U + ^ » _ g±y„ + PU y. 

x(e FC ^ ) e i;,(V, " ) ex P (Er=ia- fe (X| + a* fe y,*)ex P (-Er=i«fcn*))l 

x le^e^' ex P (Er=i«-fc(^ + c^Y*))). (4.16) 
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The Kronecker's delta 5 ee * in eq. (|4.15| ) takes a value 1 in the case J2 e i = J2 e i an d vanishes 
otherwise. Next, we compute J({d}, {w a }] {(J }, {w^}). We insert the completeness relation 



id 



nM£fcd£fc 
H W eip 



k 



E ex P E 



A; 



7 00 fc 



x|l®e A )(l(g>e" A |exp ^ 



'[2fc][fc] 



(4.17) 



between e F ^ W ' ' and e^ V and between e^ (V) and e^ (W) in the last bracket of eq. ( 4.16| ). Then 
we have the following integral expression, 



J({Q,{™ a };{C},M}) = Gx 



x exp 



x exp 



4e 



k 



d&d&d£d£ 



+ E[£*'£ fe '£fc>£fe]^ fe 
fe=i 



E 

fc=i 



[2k][k] ( 1 



v 2 _i_ a fc v* 2 



k+ 1 k ( > 



where the matrix A k and the vector are given by 
Ak= 1 " , St 



07, 


1 





-1 


1 


















1 


-1 





1 





-a^Pk- + X k + a k 1 Y k 
(3* k+ + X* + atY k * 



Using the formula for the Gaussian integral 

k 



x exp 



4e 



1 [2fc][fc] 



+ E[^'^ fe '^'^]^ fe 



fc=i 



fl (" det A)" 1/2 exp (\ £ fl* ^ , 

fc=i \ fc=i / 



in eq. ( 4.18 ), we get 

J({6>. {<».}; {<?}.{««) 



AT 

5 exp 



.jfe=l 



-(at((3 k+ + - Pt-)Y k * + a fc X fc (X,* + a* k Y k *) 

+Y k (x k + + Y k *(x* k + «j;y fc )' 



(4.18) 



(4.19) 



(4.20) 



(4.21) 
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We have used the definition of the normalization factor G (= g N /(W* • id • V*| V • id • W)). From 
eq. ( 4.15 ) with eq. ( 4.21| ) we arrive at 

Pel,...,e* N ;e N ,...,e 1 ({(j}; {Ci}) 

= ^g N nn/ c £, ^^c m ,{ Wa })c*at*},{ W * b }) 



a£Ab£B JCa JC b 



N 

x<n(-^ 2 )n-r 2 



2(Eti^ +1 ) 



X 



\j=l beB 

-5/2,+ ;^-ll/2, + V^a;^_5/2 _ VSj W-n/l,- 



' it! )^-ll/, + K)r_ 5 /2.-(Cf)^ll/, -K) 



n n 



2,7 = 1 aeA 
f>£B 



07/2,-(Ci)^l/2,-(u>a) 
^-5/2, + (C 2 )V'l/2, + K' 1 )^-5A-(C 2 )V'l/ 2 ,-(^a 1 ) 

^-i/2,-(c;~ 2 )^/2,-K" 1 ) 

03/2,+ (Cr 2 )V'5A+K- 1 )^/2,-(Cr 2 )^5/2 > -K- 1 ). 

1/2 



x (r xx ,_ 6 ({0}, {^a})ry*y., 2 ({C;}, K» 



xr X x*,-6({Ci}, H}; {Cj }, {^6}) r xr*-4({Ci}, {^a}; {(] }, M}) 
X7yy., 2 ({Ci}, W; {C;}, K})7yx*,o({Ci}, W; {C;}, K» 

xr xy , ({Ca, {^a})rx*y*,o({C;}, K»- ( 4 - 22 ) 

The functions a ±(C)i ^a±{w)i T ({(}i{ w }) are given in Appendix |C|. The integration contours 
C a and C£ encircle the following points, 

C a ■ q\firg 1 ,q 6 rg 1 -,0,±q 2 ,±q 4, ]W a >,q~~ 2 'Wa>) 

6 * 1 4 * 1 4 * * —2 * —4 * 

q w b ,q w b ,q w b ,w b ,q w b ,q w b ; 

(b e B,a' € A/{a} and i < a), (4.23) 
C 6 * : q 2 C* 2 ;r-,q 2 rZ 1 ;0,±q- 1 ,±q 3 ;w* bl ,q 2 w* bl ; 
q 6 w~ 1 ,q 4 w~ 1 ;q 4 w a , q 2 w a , w a ,q~ 4 w a ; 

(a £ A,b' £ B/{a} and j>b). (4.24) 

See eqs. ( p.6| ) and ( |C.11| ) for the definitions of txy,o arid tx*y*,o- While txy,o and tx*y*,o 
diverge at a = a' and b = b', the residues at the points q 4 Q and q 2 Q are respectively finite. The 
other functions are always finite. 

The correlation function including the local operator O is given by 

JV 

<0>=£ £ (Om) em e^ e K (4-25) 
m=l e m ,e' m =± 

E^'m = J2 Pe 1 ,...,e m ,...,e m e N ,...,e! m ,...,e 1 ({!};{!})■ (4.26) 

except for %,e^ 

For instance, the magnetization on a site m is 

<o£> = £&+-.£-- (4.27) 
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5 Discussion 



We have applied the MPA to the finite XXZ model with boundaries and calculated the 2A r -point 



function. The 2A r -point function ( 4.22 ) and the correlation function ( 4.25 ) are the main results 
of this paper. The type I vertex operator is regarded as the ZF-algebra and the boundary states 
are realized in the bosonic U q (sl2) form. The eigenstate is given by the product of them. We 
have set the vertex operators free from those physical meanings, observed those mathematical 
properties, and succeeded in analyzing the model with boundaries. Such manipulations are 
allowed only in the MPA. The vertex operators constitute the eigenstate and do not participate 
in the transfer matrix. The transfer matrix is made of the R-matrix and the boundary K-matrix. 

There are open problems. It is interesting to look for a closed equation for the correlation 
function. For the infinite XXZ model the correlation function is known to satisfy a certain 
difference equation called the quantum Knizhnik-Zamolodchikov equation |l], [7|, |26|, |27], 
and the similar equation is found also for the half-infinite case ||. In the finite chain case 
similar properties are expected. Solving that equation may give a practical method to calculate 
observable quantities. In this viewpoint it is important to investigate further the form factors 



It is also interesting to consider how to construct the excited states. In the infinite case the 
type II vertex operator plays its role. However, in the MPA, the situation is not so simple. The 
excited states are known to be obtained by the product of the ZF-algebra, the boundary states, 
and other operators that satisfy certain relations p5| . We have to find the third operators and 
show the relations among them and the type II vertex operators. 

It is possible and important to apply the MPA to other models, for instance, the massless 



XXZ model. Recently, the elliptic quantum algebra has been developed |p0| , 31] and the vertex 
operators for the degenerate elliptic algebra have been realized in the bosonic form f32]| . They 
may help us to analyze the massless XXZ model with boundaries. 

Finally, we point out a connection between the MPA and the "affinization" . In our MPA, 
the realization of the ZF-algebra has been justified by comparing with the known result. We 
like to define the ZF-algebra uniquely within a scope of the MPA. Note that the definition of 
the ZF-algebra physically means the Yang-Baxter equation for the transfer matrix on the half- 
infinite lattice. Such a transfer matrix is obtained by extending the original one into the infinite 
chain limit, i.e. by the "affinization" of symmetry of the model. We expect that the ZF-algebra 
is uniquely realized by the affinization of the transfer matrix. The MPA may be mathematically 
defined by using the idea of the affinization. 
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A R-Matrix and Boundary K-Matrix 

We summarize some basic formulas for the R-matrix and the boundary K-matrix of the XXZ 
model. 
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A.l R-Matrix 

The R-matrix is 



R(0 



<0 



T7T 



1 - qX 



(i - g 2 )C 
(i - C 2 )g 



(A.l) 



where 



«(0 



(g 4 C 2 ;g 4 )oo(g 2 r 2 ;g 4 ) 
(g 4 C- 2 ;<? 4 )oo(g 2 C 2 ;g 4 ) 



-C, (*;p) 



oo 

x - n a 

n=0 



(A.2) 



The R-matrix satisfies the Yang-Baxter equation, the unitarity and the crossing symmetry, 



i?12(Cl/C2)^13(Cl/C3) J R23(C2/C3) = 

#12(Cl/C2)#2l(C2/Cl) = 
e' e' 

A.2 Boundary K-Matrix 

In this paper, we choose the diagonal K-matrix, 

K((;r) 



^23(C 2 /C3)^13(Cl/C3)^12(Cl/C2) 

1, 



jqSHrta/Ci). 



(A.3) 
(A.4) 

(A.5) 



where 



/(C;0 



y>(CV) 



¥>(z;r) 



1-rC 2 



¥>(CV) ' ^ ' ; (9 2 r^ 4 )oc(<Z 8 * 2 ;<Z 8 )oo' 
In the text, we denote K((;r±) by K ± (C)- The K-matrix has the following properties, 

K(C;r)K(cV) = 1, 

K b a {-q^C l ;r) = ^i?^(-gC 2 )<(C;r). 



(A.6) 



(A.7) 



(A.8) 
(A.9) 



a'V 



B Vertex Operators 

We formulate the type I vertex operators. (The type II vertex operators are omitted since they 
are not used in this paper.) We consider the bosonic Fock space 

W« = C[a_i, . . .] ® (e„ eZ Ce A+ria+<5 " a / 2 ), (B.l) 

for i = 0, 1. The commutation relations among the bosonic operators {a n } (n eZ/{0}) are 



["mi O n ] — <^m+n,0 



[2m] [m] q n — q n 

m q — q 1 



(B.2) 



The operators e a and z act as 



e Q . = e a+/3 , z 9 ■ e a = z^e a ■ z d , 



(B.3) 
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where [d, a] = 2 and [d, A] = 0. These operators act also on the dual Fock space TC*^ as 

e P . e « = 6/3+°, e a • z 9 = z 9 ■ e a z^ . (B.4) 

We further define that (e a ) = and (z 9 ) = 1. 

The vertex operators are defined as the intertwiners that have the following map, 

* (1_M) (C) : -^w^gv; 

$(i-M)(f) = J2 ^- l ' l) (C) ® « e) (B.5) 

e=± 

and are realized in the following bosonic form, 

$(_!-*.*) ( C ) = e P(O e Q(O 0e «/2(_,3 C 2 ) (» + i)/2 c -< > (B.6) 

■ / dw (l-g 2 )wC 



c 2vrig(^-g 2 C 2 )(w-g 4 C 2 ) 
x e P(0+%) e Q(0+sH 8e -«/2 ( _,3 (2) O + i)/2 r a r , (B?) 



where 



p (0 = E jgp 7fc/2 c 2fc , Q(0 = - E ^~ 5k/2 <~ 2k > (b.8) 

oo oo 

= - E *M = E Ss* 72 ""* ( B - 9 ) 

fc=i M fe=i M 

The integration contour C in eq. (|B.7| ) encircles the point g 4 £ 2 but not the point q 2 C 2 - 
Similarly the dual vertex operators are defined as 

$*(i-M)(£) = £ g, v * ; (b. 10 ) 

where v* is the dual vector of v that satisfies v* ■ v e > = 8 ee i. The dual vertex operator is related 
with the original one through 

$*(l-M)(£) = ^^(-q^Q. (B.ll) 

The label (1 — i, i) of the vertex operator is omitted if there is no fear of confusion. The following 
relations are known, 

E ^(Cl/C2)* ei (Cl)*e 2 (C2) = *4(C2)^ 1 (Cl), (B.12) 
ei,e 2 =± 

^$:(c)^(o = id (B.i3) 

where 5 is a scalar constant, 

5 = (R14) 
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Formulas 



summarize the formulas that are used for the derivation of eq. ( 4.22| ) from eq. 



,±(z) = exp [f^ [ ^-q ak z k p k± 



' {q 2a ~ l z 2 - g 8 )oo(g 2a - 3 ^ 2 ; 9 8 )oo \ V2 {q a+1 ' 2 r ± z- g 4 



Ik' 



\fc=i 

l _ q 2a-7 w 2\ V2 , 



1 - q 2a ~ 5 w 2 I 1 - g a - 5 /2 r± u; ' 



' (^+3^2. g 8 )oo(g2a+ 9 z 2. g 8 )oo \ V3 (g «+5/2 r±z . g 4 



((? 2 a+ 7 z 2. g 8 )oo(g 2 a+ 5 z 2. g 8 )oo J ( q a+l/2 r±z . g 4) 



oo 



C±H = exp ( ]T 



,2o!+l,,,2\ 1 / 2 



1 - q^ a+L w^ \ ' 1 



X _ q 2a+3 w 2 ) x _ „a-l/2 



rxxAidhM) = exp ■ ? afe J fe I,J 

n (9 a+10 C 2 C|;g 4 )oo(i - g a+4 CNaO(i - g a+4 Qa) 
11 11 (^+8c 2 C^;g 4 )oo(l-g Q WaO(l-^ +2 ^aO ' 



i,i'=l a,a'GA 



= tt n (g" +5 C 2 C7 2 ; g 4 )oo(i - g^dVXl - rV%) 

itoitl (^ +3 C 2 C,7 2 ;9 4 )oo(l - q-w a w-, l )(l - q 
Txx'Aiti}* M {(*}, {<}) = exp ( f; 

A 7 " 

= n n 



(g Q+8 C 2 C- W)oo(i - g" +4 C 2 <)(i - g a+2 c; 2 ^a) 

,, ,„_ A {q a+& C,Kf\Q A Uil-q a WaWt)(l-q-+ 2 w a wt) ' 

rxY^aa, W;{c;},K>) = exp (f; 



V.k=l 



A? 



=! M (^ +4 C 2 C*- 2 ; g 4 )oo(i - gVT^ - f+vr 1 ) ' 
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TV*.,«({Ci}, W; {C*}, K» = exp (f) [ ^q ak Y k Xt 



\k=l 



N 



r^,a(te},{^};{c;},K}) = exp 



7X.y.,a({C?}, K» = exp ( £ 



\fc=i 



TY*y*A{Ch {<}) = exp ( E M^y, 



\fc=i 



JV 



n n (g" +2 Cr 2 C;r 2 ; g 4 )oo(l - g"Cr 2 <' 1 )( 1 ~ g a Cf 2 <' 1 ) 
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